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PACS 73.20.Mf - Collective excitations (including excitons, polarons, plasmons and other charge- 
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Abstract. - Surface plasmon polaritons (SPPs) have recently been recognized as an important fu- 
ture technique for microelectronics. Such SPPs have been studied using classical theory. However, 
current state-of-the-art experiments are rapidly approaching nanoscales, and quantum effects can 
then become important. Here we study the properties of quantum SPPs at the interface between 
an electron quantum plasma and a dielectric material. It is shown that the effect of quantum 
broadening of the transition layer is most important. In particular, the damping of SPPs does 
not vanish even in the absence of collisional dissipation, thus posing a fundamental size limit for 
plasmonic devices. Consequences and applications of our results are pointed out. 
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The excitation and propagation of surface modes at 
plasma interfaces have long been important in space 
physics, magnetic confinement fusion, and laboratory 
plasma physics [1-4]. Moreover, studies of electron os- 
cillation excitations (surface plasmons, surface plasmon 
polaritons (SPPs) and magnetoplasmons) [5-7] in nano- 
structured systems [8-11] have recently attracted much 
interest. It has been found that at condensed matter inter- 
faces, such plasmon excitations could be of crucial impor- 
tance in future electronic components [13-16], the latter 
■ referred to as plasmonic devices [17, 18]. 

Furthermore, the field of quantum plasmas has recently 
developed rapidly [19]. This field started already in the 
1960's, when Pines studied the excitation spectrum of 
quantum plasmas [20] , with a high density and a low tem- 
perature as compared to normal plasmas. In such sys- 
tems, the finite width of the electron wave function makes 
quantum tunnelling effects crucial, leading to an altered 
dispersion relation. Since then, a number of theoretical 
studies of quantum statistical properties of plasmas have 
been published (see, e.g., Ref. [21] and references therein). 
For example, Bezzerides & DuBois presented a kinetic the- 
ory for the quantum electrodynamical properties of non- 
thermal plasmas [22], while Hakim & Heyvaerts used a co- 
variant Wigner function approach for relativistic quantum 
plasmas [23]. It has also been shown that, under certain 
conditions, plasmas can display remarkable properties due 
to the quantum properties of the constituents. Thus, there 
are quantum multistream instabilities [24,25], quantum 



modified Zakharov dynamics [26, 27] together with soli- 
ton formation and nonlinear quantum interactions [28,29], 
spin effects on the plasma dispersion [30,31], quantum 
plasma turbulence [32], ferromagnetic plasma behaviour 
and Jeans-like instabilities due to quantum effects [33]. 

Many of the current studies involving quantum plasmas 
are motivated by the rapid experimental progress and de- 
velopment of new materials, e.g., nanostructured materi- 
als [34] and quantum wells [35], and the laboratory real- 
ization of ultracold plasmas [36] and experimental demon- 
stration of quantum plasma oscillations in Rydberg sys- 
tems [37]. Quantum dispersive effects can also be impor- 
tant for diagnostics of inertial fusion plasmas [38] . In par- 
allel, the field of plasmonics and its use of surface waves, 
such as surface plasmon polaritons (SPPs), has emerged 
as a new route to electronic devices [12]. 

In this Letter, methods from quantum plasma physics 
are used for analyzing SPPs in nanoscale systems. In par- 
ticular, we determine the dispersion relation for quantum 
SPPs on a conductor-dielectric interface. It is shown that 
wave function dispersion introduces an intrinsic damping, 
even in the absence of collisions. Such damping is due to 
the irreversible propagation of resonant plasmons towards 
lower density regions. This is of importance for the short 
length scales in forthcoming electronic components based 
on SPPs in waveguide slots. An expression for the damp- 
ing length, limiting the size of such devices, is presented. 
For wavelengths of the order of tens of nanometers, the 
propagation distance for the SPPs is only a fraction of the 
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wavelength. 

The governing equations for the electrostatic electron 
dynamics are [19] the continuity equation dtn + V ■ (nv) = 
0, the momentum conservation equation 



mn(9t+v-V)v = -enE-Vp-mnw- 



h 2 n^ 
2m 



(1) 

and Poisson's equation V • (e E + P) = Jij — n, where n 
is the electron density, m is the ion density (here we treat 
the lattice constituents in terms of the ion density, in or- 
der to demonstrate the analogy to classical plasmas and 
surface plasmons), m is the electron mass, v is the elec- 
tron velocity, e is the magnitude of the electron charge, E 
denotes the electric field strength, p is the electron pres- 
sure, h is Planck's constant, eo is the vacuum dielectric 
constant, v is the collisional frequency between the elec- 
trons and the ions, P = eo(ed — 1)E is the polarization 
due to bound charges in the system, and ed is the per- 
mittivity. The last term in Eq. (1), the gradient of the 
Bohm-de Broglie potential, corresponds to the effect of 
wave function dispersion. 

We consider the ions as stationary, and write the per- 
turbed electric field according to <5E = — V56. Moreover, 
we assume that the conducting medium is semi-infinite 
in the positive x-direction, so that the stationary back- 
ground electron density no is a function of x (see Fig. 
1). Next, we linearize the governing equations. In the 
cold classical case we obtain the SPP dispersion relation 
V • [e(x)V86] = 0, where e(x) = Cd(x) — uj 2 (x)/uj(uj + iv) 
is the total dielectric function for a cold classical plasma 
with collision frequency v, and u> p (x) = [e 2 no(x) /eom] 1 / 2 
is the electron plasma frequency. We note that the colli- 
sional frequency is a function of the ion density. In the 
general case, we obtain 
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where we use the index notation V — ► d a etc., as well 
as Einstein's summation convention and let all perturbed 
quantities be of the form g(x) exp(iky — iut) [39]. Here 
the dielectric tensor operator is given by 
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speed. 
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where c s = [(dp/dn)\, 
We will from here on consider 
(47r 2 ?i 2 /5m)(3/87r) 2 / 3 7j 5 / 3 of a non-relativistic degener- 
ate electron gas, so that the sound speed becomes c s = 
(27rVV3m)(3n /87r) 1 / 3 . 

We consider a plasma which can be divided into three 
regions (see Fig. 1). For x < 0, the unperturbed plasma 
density is zero. For < x < d, we have a transition layer 
where the plasma density is a monotonously increasing 
function of x. Lastly, the plasma density is constant for 
x > d. Furthermore, for the surface waves to be only 
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Fig. 1: The structure of the electron density profile for the SPP 
excitation. The width d of the transition layer is determined by 
quantum effects. A dielectric material, with electron density 
rid(x), is on the left, while the metal with the corresponding 
plasma oscillations is on the right. The small panel in the 
upper left corner shows the Airy-like structure of the density 
perturbations Sn/no- 

weakly damped, the condition kd ^ 1 must be fulfilled. 
The physics inside the transition layer is then essentially 
one-dimensional, since k -C d/dx applies here. 

Integrating Eq. (2) across the transition layer, we obtain 



c d x 86 + ike xv i 
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Moreover, 86(d) and 86(0) are related by 



(4) 
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dx 
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no(x)v x (x) dx (5) 



The integrals in Eqs. (4) and (5) are correction terms pro- 
portional to kd. Similarly the off-diagonal quantum terms 
in the left hand side of (4) are also proportional to one 
of two small parameters, namely k^f^/m^Lu 2 (from the 
Bohm-de Broglie potential) or k 2 c 2 /uj 2 (from the Fermi 
pressure). If we drop both these quantum terms and let 
d — > 0, we obtain the wellknown electrostatic dispersion 
relation to = Wp (1 + ^)^ 1 ^ 2 , where — e<i(x < 0) 
is the constant dielectric permittivity for x < 0, and oj p °^ 
is the plasma freqency at x = d. Within this approxima- 
tion, the electrostatic surface wave has a zero damping as 
well as a zero group velocity. The quantum terms modify 
this result, since d now remains finite. To adress this case 
properly, we must determine no(x)v x (x) inside the transi- 
tion layer. Using the approximation k -C d/dx inside this 
region we find 



Ma;K(aO] = ^(0)^(0). 



(6) 



We can here treat 86(x) as a constant within the 
layer, in which case the right hand side of (4) reduces to 
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5<j)(0) J Q k 2 e(x) dx. Moreover, when the inhomogeneous 
derivatives on the left hand side of (4) vanish, the remain- 
ing quantum terms become proportional to V 2 8<t>. How- 
ever, for k 2 c 2 /oj 2 <C 1 and k 4 H 2 /m 2 uj 2 <C 1, the homoge- 
neous regions are classical to first order, i.e., V 2 <5</> = 0. 
The dispersion relation for SPPs then becomes 



,(°) 



(1 + 4V 



k 

1 + 2 / ^ + q ^ dx 
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where q(x) is the solution to Eq. (6) for the plas- 
mon field inside the transition layer, defined by q(x) = 
ieno(x)v x (x)/u)kS(p(0). 

The integral contribution in (7) depends on the equi- 
librium profile. Before turning our attention to the full 
quantum equilibrium profile we consider the following. 

i) Classical case. Here a finite width of the elec- 
tron distribution in the transitions layer is determined 
by a finite width of the ion profile. From (6) we ob- 
tain q[x) = e(0)/e(x). The singularity in this equa- 
tion can be dealt with using the Landau prescription, 
which results in an imaginary contribution to the integral, 
Im{/ qdx} = ine(0)(de/ dx)~2 in Eq. (7), where the sub- 
script indicates that the derivative should be evaluated at 
the resonant surface where e = 0. The imaginary part 
representing damping of the SPP due to energy transfer 
to the resonant plasmons in the transition layer. 

ii) Fermi pressure effects. Assuming that the width is 
determined by the ion-density profile we drop the Bohm- 
de Broglie potential. In this case, (6) has the solution 
q(x) = Gi(xde/dx) + iAi(xdefdx), where we for simplic- 
ity assume that e(x) is a linear function. Here Gi and Ai 
denote solutions to the inhomogeneous and homogeneous 
Airy equation [40], respectively. The proper causual solu- 
tion has been found by taking the nondivergent solution 
to Eq. (6) that represents resonant plasmons propagat- 
ing towards lower densities. Due to negligible reflection 
of the plasmons at the dielectric interface, the energy is 
irreversibly lost, even in the absence of dissipation. As 
a consequence, the damping of the surface wavecoincides 
between cases i and ii. 

Next, we investigate the cquibrium density profile in the 
transition layer, determined by 
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where we have introduced the normalized distance 
x = x/^h/V^mujp^) 1 / 2 , the normalized density n = 
n (x)/n§\ taken the normalized ion-density as a step 
function H(—x), and introduced the dimcnsionless pa- 
rameter n = 47r 2 (a/V2m^ 0) )(3n[ ) 0) /87r) 2 ' /3 . The relative 
importance of the the Fermi pressure and the Bohm dc 
Broglie potential is determined by the parameter rj. We 
first consider the case i] <C 1, i.e., the Bohm-de Broglie 
potential dominates over the Fermi pressure. Noting that 
n — H(—x) is of order unity after passing the ion density 



step, the normalized width d must be of order unity, i.e., 
in dimensional units we have d ~ (h/mco^) 1 / 2 . On the 
other hand, assuming a dominating Fermi pressure, i.e., 
n » 1, the width becomes d ~ (jjfi/mwp"') 1 / 2 . For a solid 
state plasma, which is the case of most interest, r\ is of or- 
der unity, and thus we get d ~ (Ti/moj^) independently 
of which of the two estimates we use. Assuming 77 <~ 1, we 
can deduce a number of properties of the resonant plas- 
mon field inside the transition layer using the solutions to 

(6) . 

With no (a;) determined by (8), we can now treat the 
full quantum case. We start with the simple observation 
that q(x) is independent of k, as terms involving k can 
be neglected in (6) due to k <C d/dx. Solving Eq. (6) 
to find a full solution for q(x) is a nontrivial task, requir- 
ing a solution to the equilibrium equation (8), that has to 
be found numerically. However, a qualitative understand- 
ing can be obtained by Taylor expanding Eq. (6) close to 
the classical plasmon resonance e(x) = 0. We obtain an 
Airy like behavior of the plasma oscillations along the re- 
direction in the transition layer, depicted in the upper left 
panel of Fig. 1. The full quantum case gives two major 
changes to the classical and Fermi pressure cases. Firstly, 
the inhomogeneitics in the transition layer shifts the plas- 
mon resonance somewhat towards higher densities. Sec- 
ondly, the typical plasmon wavelength in the x-direction 
becomes comparable to the width of the transition layer. 
Using this as a basis for an estimate, together with the 
condition that no plasmons propagate towards higher den- 
sities in the transition layer, we find that f Q d q(x) dx w 
(-Kh/muj) 1 / 2 (0.6 + 2i). The uncertainty in this expression 
is directly linked to the uncertainty in the unperturbed 
density profile and the expression for the transition layer 
width d. Making a similar estimate of the term J Q e(x) dx 
in (7), we use a linear profile for e(x), in which case we 
obtain t(x) dx — (-Kfi/muj) l / 2 [t(x = d) — e(x = 0)] w 0, 
where the last equality follows from the dispersion relation 

(7) . A more accurate density profile gives a finite contri- 
bution. However, the plasmon coupling term <x q(x) dx 
will still dominate. Thus, we can write our final dispersion 
relation for the SPPs as 
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The group velocity and the damping rate of the electro- 
static SPPs are V g — duo/dk = 0.6V q and Imw = 2kV c 



respectively, where we have introduced the characteristic 
velocity V q = [Tiuj^ /m(l + e^)] 1 ^ 2 - In Fig. 2 we display 
the real part of the dispersion relation (9). Apart from the 
shortest wavelengths, the dispersive properties are domi- 
nated by electromagnetic effects. In the short wavelength 
limit the group velocity of the classical dispersion relation 
approaches zero, whereas the quantum corrected value V g 
approaches a small but nonzero constant value V q . Let 
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Fig. 2: The classical dispersion relation (thin curve), and the 
quantum dispersion relation (thick curve) for the SPPs. The 
dashed horizontal line gives the classical resonance frequency. 
For wavenumbers larger than k c = (e^Wp ^) 5 ^ 6 (m/7ic 4 ) 1 / 6 , the 
quantum contribution to the group velocity is larger than the 
classical electromagnetic contribution. We note that quantum 
damping, as represented by Imui in (9), can be important for 
wavenumber lower than k c . 

us point out that there might be classical effects not in- 
cluded in our model that can modify the given picture, e.g. 
thermal effects which can induce a fc-dependence of e. We 
note, however, that within a fluid model, the density per- 
turbation of the surface wave approaches zero, which leads 
to negligible modifications of our model. Within a kinetic 
picture, on the other hand, thermal effects induce classi- 
cal modifications with a fc-dependence of e, which might 
lead to modifications of the dispersive properties also in 
the short wavelength limit. 

However, the main importance of the dispersion relation 
(9) concerns the dissipative (i.e. imaginary) part, which 
can have significant consequences for the rapidly emerg- 
ing fields of plasmonics [15,17,18], where plasmon wave 
propagation along metallic-dielectric interfaces is studied 
as a means to pave the way for even smaller and faster 
electronc circuits [15,16]. Naturally, energy losses must 
be minimized if this undertaking is to be successful. Cur- 
rently the propagation distances of the SPPs range from a 
few 100 nanometers up to tens of microns [15]. The way to 
improve performance has so far been to consider wavguide 
slots [15, 17, 18], where most of the energy is distributed in 
a dielectric, limiting the collisional losses in the metal sur- 
face. By cooling the system collisional losses can be further 
minimized [41]. However, when reducing the size towards 
the nanoscale, the wavelength of the SPP is also decreased. 
Assuming that collisional effects are minimized when ap- 
proaching such a regime, our treatment of quantum losses 
becomes important. For such nano-sized system, a damp- 
ing of the order Im(w)/|w| ~ k(?i/muj v ) 1 l' 2 , as given by 
Eq. (9), becomes crucial, and sets a fundamental limiting 
factor for how small systems that can be designed. The 
damping length of the SPP due to quantum effects is given 
by <5 S p = V^cm/Imw, where V g , em « (e { ^ / c) 2 {uj^ / k) 3 is 
the group velocity including electromagnetic effects [4] in 
the short wavelength region, while Imw is given by (9). 



It should be stressed that the electromagnetic contribu- 
tion to the group velocity dominates over the quantum 
induced contribution V q for wavelengths A > 30 nm. As- 
suming that the dielectric consists of Si02, we have [42] 
e d°^ ^3 — 5, and with the plasma frequency of the metal 
as u>p «4x 10 15 s _1 , we obtain 

\ 100 nm J 

where A = 2n/k is the wavelength. Thus, due to the strong 
wavelength dependence in (10), /im- waves can propagate 
without significant quantum damping, while decreasing 
the scale much below the /im regime will affect the ef- 
fective propagation distance. For example, for A ~ 30 nm 
the damping length #sp ~ 10 nm. Although different ge- 
ometries may affect the possibilities to design smaller de- 
vices [15], our result (10) is robust, and its consequences 
must therefore be considered in the design of plasmonic 
devices. 
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